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Abstract 

Organisms with coexisting sexual and asexual populations are ideal models for 
studying the consequences of either reproductive mode on the quantitative genetic 
architecture of life-history traits. In the aphid Rhopalosiphum padi, lineages differ- 
ing in their sex investment coexist but all share a common parthenogenetic phase. 
Here, we studied multiple genotypes of R. padi specialized either for sexual and 
asexual reproduction and compared their genetic variation in fitness during the 
parthenogenetic phase. Specifically, we estimated maintenance costs as standard 
metabolic rate (SMR), together with fitness (measured as the intrinsic rate of in- 
crease and the net reproductive rate). We found that genetic variation (in terms 
of broad-sense heritability) in fitness was higher in asexual genotypes compared 
with sexual genotypes. Also, we found that asexual genotypes exhibited several 
positive genetic correlations indicating that body mass, whole-animal SMR, and 
apterous individuals production are contributing to fitness. Hence, it appears that in 
asexual genotypes, energy is fully allocated to maximize the production of partheno- 
genetic individuals, the simplest possible form of aphid repertoire of life-histories 
strategies. 
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Introduction 

One of the most important attempts to formalize a general 
theory in biology was Fisher's fundamental theorem of nat- 
ural selection (FTNS), according to which the instantaneous 
rate of increase in the mean fitness of a population equals 
its additive genetic variance in fitness (Fisher 1930; Okasha 
2008; Bijma 2010). However, the realism of the theorem has 
been questioned as the FTNS predicts a continuous increase 

1414 



in fitness maintaining all other things equal (including the 
environment and the average effects of alleles) (Bijma 2010). 
Among other things, nonadditive genetic variation and re- 
source limiting leading to competition and deterioration of 
the environment are common factors in natural populations 
that violate these assumptions. Then, the average fitness of the 
population cannot always increase (Okasha 2008). The initial 
and perhaps most popular empirical support for the Fisher 
fundamental theorem was based on the fact that life histories 
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have low heritability because they are constant targets of 
directional selection (Gustafsson 1986; Roff and Mousseau 
1986; Mousseau and Roff 1987). However, this was later de- 
bated because of the neglected contribution of environmen- 
tal variation that life-history traits accumulate (Price and 
Schluter 1991; Merila and Sheldon 2000; Kruuk et al. 2003). 

The modern interpretation of the FTNS, which considers 
the partial change in mean fitness under a constant environ- 
ment (Okasha 2008), has renewed the interest in the quantita- 
tive genetics of fitness traits (Reale and Festa-Bianchet 2000; 
Crow 2002; Bolnick and Nosil 2007; Okasha 2008; Ewens 
2011; Guillaume 2011; Hendry et al. 2011). The FTNS can 
thus be used as a null hypothesis for the empirical study of 
adaptation, permitting to test its assumptions in experimen- 
tal settings (Reale and Festa-Bianchet 2000). With this inter- 
pretation, empirical supports for the FTNS became varied 
(Kruuk et al. 2000; Reale and Festa-Bianchet 2000; McCleery 
et al. 2004; Bolnick and Nosil 2007). 

Although several authors already advocated the use of 
cyclic parthenogenetic organisms for the study of the fitness 
consequences of sexual reproduction (Colegrave et al. 2002), 
few experimental studies have addressed in these organisms 
the problem of genetic variation in fitness, experimentally 
(Pfrender and Lynch 2000; Morgan et al. 2001). In cyclically 
reproducing organisms, sexuality is a phenotypically plas- 
tic trait that can be triggered by environmental conditions 
(De Meester et al. 2004). This is the case of aphids, for which 
the alternance of parthenogenetic (asexual) and sexual repro- 
duction has important consequences on the genetic structure 
and architecture of natural populations (Simon et al. 1991; 
Simon et al. 2002; Halkett et al. 2008; Nespolo et al. 2008; 
Nespolo et al. 2009). In fact, parthenogenesis is the main 
reproductive mode of aphids while sexual reproduction is 
triggered once a year to anticipate stressful conditions of 
cold winter. During warm periods aphids reproduce asexu- 
ally (producing viviparous clonal offspring) exhibiting high 
population growth rates (Rispe et al. 1996). When the cold 
season begins, sexual morphs are produced, which mate and 
lay frost-resistant eggs (Rispe et al. 1996). Thus, the main 
known advantage of sexual reproduction is the resistance to 
freezing conditions, which typically kill the asexual forms 
(Newton and Dixon 1988; Simon et al. 1991; Rispe et al. 
1996; Simon et al. 2002; Nespolo et al. 2009). 

In its distribution range, the bird cherry-oat aphid, 
Rhopalosiphum padi (L.) displays genotype specialization 
associated to reproductive mode variation (Simon et al. 
1991; Hulle et al. 1999), which has consequences on its life- 
history patterns, dispersal traits, and reproductive investment 
(Simon et al. 1991; Simon et al. 2002). Two kinds of geno- 
types can be identified in terms of their capacity for producing 
sexual morphs: "asexual" (i.e., sustained asexuality without 
a full commitment to the sexual reproduction) and "sexual" 
(i.e., obligate and regular alternation of asexual and sexual re- 



Reproductive Specialization Costs and Genetic Variance 

production) genotypes. Sexual and asexual genotypes share 
a common summer parthenogenetic phase on cereals and 
other grasses. However, they live on separate hosts the re- 
maining time because sexual reproduction is associated with 
host shift in this aphid. Earlier studies showed that asex- 
ual lineages may originate from sexual ones through diverse 
mechanisms so that they form a mixed array of phylogeneti- 
cally distinct genotypes (Delmotte et al. 2003). Hence, -R. padi 
constitutes a good system to study the impact of reproduc- 
tive mode variation on the genetic architecture of life-history 
traits because it allows within and between comparisons of 
sexual and asexual lineages with distinct genetic backgrounds 
(Delmotte et al. 2003; Nespolo et al. 2008). While the popula- 
tion genetics and phylogeny of sexual and asexual genotypes 
of R. padi have been addressed in several studies, little at- 
tention has been paid to their respective ecology and fitness. 
In particular only few works have studied the life histories 
of reproductive specialization in R. padi. In previous works, 
Nespolo and co-workers determined the effects of tempera- 
ture on asexual and sexual morph production (Nespolo et al. 
2008) and the life-history consequences of phase induction 
(Nespolo et al. 2009). Both studies were performed on the 
basis of genetic variances in several traits (i.e., comparing 
"G matrices"), and showed that sexuality and reproductive 
mode variation have profound effects on the structure of 
the G matrices. However, both studies analyzed an undif- 
ferentiated array of sexually and asexually specialized geno- 
types. Given the strong and consistent pattern of reproduc- 
tive specialization that _R. padi exhibits (Halkett et al. 2006), 
it is expected that during the shared parthenogenetic phase, 
asexually specialized genotypes would be benefited compared 
with sexually specialized ones, since the latters are benefited 
in winter (Moran 1992; Rispe et al. 1996). In fact, previous 
studies have indicated that asexual genotypes maximize their 
survival and fecundity under these conditions, which can be 
interpreted as being "adapted" to such (summer) conditions 
(Rispe etal. 1996; Hulle etal. 1999). The FTNS tells "the rate of 
increase in fitness of any organism at any time is equal to its ge- 
netic variance in fitness at that time" (Edwards 1994). Hence, 
if asexual genotypes are in the process of being "adapted" 
to summer conditions, its average fitness should be increas- 
ing (i.e., genetic variance or heritability of fitness should be 
higher than zero). On the other hand, if asexual genotypes are 
already adapted, being fixed all their alleles, then average pop- 
ulation fitness should not longer increase (i.e., genetic vari- 
ance or heritability of fitness should not different than zero). 
Seasonality and population subdivision would promote the 
former conditions, as populations would never attain the re- 
quired equilibrium conditions for the adaptation to be fixed. 
During summer, sexual genotypes would be at disadvantage 
compared with asexual genotypes, and hence the average 
increase of population fitness would be very low or unde- 
tectable (i.e., genetic variance or heritability of fitness should 



©2012 The Authors. Published by Blackwell Publishing Ltd. 



1415 



Reproductive Specialization Costs and Genetic Variance 



M. J. Carter ef al. 



not different than zero). Also, several specific features of the 
adaptation process would appear in the asexual genotypes, 
compared with the sexual ones. For instance, the allocation of 
energy to less costly asexual phenotypes could be translated 
into lower energy costs. However, it is not clear what to ex- 
pect during the asexual phase, where sexual genotypes are not 
constrained to produce sexual forms. One rationale would 
be that the genetic program of sexual genotypes maintains 
its constraints even when not producing sexual forms, thus 
reproducing at a higher energy cost during the partheno- 
genetic phase, compared with asexual genotypes. Hence, 
higher energy metabolism would be expected in sexual 
genotypes. 

It is important to note that differential patterns of energy 
allocation, could be reflected either in direct energy measure- 
ments (i.e., metabolic rate) or in differences of body mass (i.e., 
somatic tissue). In order to test these predictions, we com- 
pared asexual and sexual genotypes during the asexual phase, 
in _R. padi. 

Material and Methods 
Biological system 

Our biological model was the bird cherry-oat aphid 
(R. padi; Fig. 1). As typical for aphids, "sexual" genotypes of 
_R. padi reproduce through cyclical parthenogenesis, whereby 
a parthenogenetic phase on many species of Poaceae (cereals 
and other grasses) alternates with a sexual phase on a specific 
winter host, which is the tree species Prunus padus L. (Simon 
et al. 1991; Rispe et al. 1999). The shift from parthenogenetic 
to sexual reproduction occurs by the end of summer and is 
triggered by increasing night length and decreasing temper- 
ature (Dixon and Glen 1971). By contrast, asexual genotypes 
reproduce all year long parthenogenetically on a wide range 
of Poaceae plants. Sexual and asexual genotypes generally co- 
exist in regions where cold and warm winters fluctuate, which 




Figure 1. Model organism Rhopalosiphum padi during parthenogenetic 
(viviparous) reproduction. Photo by Bernard Chaubet. 



is the case of western France (Simon et al. 1991; Halkett et al. 
2008). 

From end of fall to late spring, asexual and sexual geno- 
types are clearly separated in space, being on Poaceae and 
on P. padus, respectively. This spatial separation conveniently 
allows ascertaining sexual and asexual genotypes since no 
morphological characters exist to distinguish between them. 
Asexual genotypes were sampled by end of fall beginning of 
winter on barley or wheat plants separated by at least 1 m in 
order to minimize collection of aphids from the same clonal 
genotype. Sexual genotypes were collected on P. padus af- 
ter they had just hatched from eggs, ensuring that they did 
not undergo asexual reproduction. To initiate clonal (isofe- 
male) genotypes in the laboratory, parthenogenetic females of 
asexual and sexual genotypes of R. padi were collected respec- 
tively in Pleine-Fougeres (Brittany, France) in the vicinity of 
Mont-Saint Michel in autumn 2006 and 2007, and in La 
Chapelle Thouarault (Brittany, France) in the vicinity of 
Rennes in spring 2008. Each parthenogenetic female was kept 
isolated on a wheat plant and allowed to establish an isofemale 
genotype. Asexual and sexual genotypes reared in the labora- 
tory were then genotyped at eight microsatellite markers to 
characterize their multilocus genotype or MLG (see Halkett 
et al. 2008 for methodological aspects). Fifteen asexual geno- 
types and 20 sexual genotypes were then selected for the 
experiments based on their MLG (Table Al). Because clonal 
copies of the same genotype are frequent in asexual lineages of 
R. padi (Delmotte et al. 2002), we were constrained in the 
number of asexual genotypes. Genetic analysis also permit- 
ted to check that aphid genotypes used here did not belong 
to differentiated genetic clusters associated with either repro- 
ductive modes but rather represented an array of sexual and 
asexual lineages with distinct genetic backgrounds (Fig. Al). 
Each genotype was reared in isolation under conditions en- 
suring continuous parthenogenetic reproduction, during ca. 
20 generations (16 h light:8 h dark, 18°C) prior experiments. 
Wheat was used for both rearing and experiments. 

Fitness measurements 

Sexual and asexual genotypes were monitored in controlled 
conditions to assess various life-history traits. Experimental 
conditions for fitness measurements were set at 18°C and 16 h 
of light to ensure sustained parthenogenetic reproduction in 
both sexual and asexual genotypes of .R. padi (Dixon and 
Glen 1971). To isolate maternal and grand-maternal effects, 
we started rearing parthenogenetic females of each geno- 
type in isolation, two generations before the experiments. 
For trait measurements, we kept the clonal genotypes at the 
rearing conditions and we assessed, on eight replicates per 
genotype, the following traits: age at maturity of focal indi- 
viduals (number of days elapsed from their birth to the first 
offspring produced) and fecundity of winged (=production 



1416 



© 2012 The Authors. Published by Blackwell Publishing Ltd. 



M.J. Carter ef al. 



Reproductive Specialization Costs and Genetic Variance 



of winged individuals from the focal individual) and wing- 
less (=production of wingless individuals from the focal 
individual). With this information, we constructed life ta- 
bles and estimated fitness as the intrinsic growth rate (r), by 
means of the discrete version of the Euler-Lotka equation, 
Yl"=w e_rl ii m i = 1> and by computing lifetime reproduction 
of each focal individual, as the total number of individuals 
produced during its lifetime. 

Standard metabolic rate (SMR) 

As indicated for animals with very low rate of metabolism, 
rates of carbon dioxide production (VCO2) were used as 
measure of SMR (Lighton 2008). VCO2 was determined us- 
ing "closed system" metabolic chambers (Lighton 2008), con- 
sisting of glass 2.5 mL hermetic syringes fitted with three-way 
valves (see also Chappell 1983; Ashby 1997). Animals were 
weighed (adult body mass = Mb) to the nearest micrograms 
with a microbalance and then placed, individually, in the 
syringes. The syringes were sealed from the atmosphere and 
placed in a temperature controlled, dark incubator for the du- 
ration of the measurement period (ca. 90 min). Three blank 
syringes served as controls for each series of measurements. 
We injected the air of the syringe into a tube connected to the 
CO2 analyzer after passing through C0 2 -absorbent granules 
of Baralyme™ and Drierite™. At the end of the measure- 
ment interval, CO2 concentrations were determined using 
an infrared CO2 analyzer Li-Cor 6251 (Lincoln, NV, USA) 
capable of resolving difference of 1 part per million (ppm) of 
CO2 in air. We computed the total amount of C0 2 produced 
per aphid as the area below the curve of CO2 concentration 
in time. Then, the measures were transformed from ppm to 
mL CO2 h _1 , taking into account the flow rate and incuba- 
tion time for each syringe. Therefore, each measurement is 
an average of CO2 production over several hours. Accord- 
ing to previous measurements, technical errors associated 
with this measurement method are small (Castaneda et al. 
2010a,b), and its simplicity allows simultaneous measure- 
ments of a large number of individuals. In this experiment, 
we measured SMR for a total of ca. 312 individual aphids. 

Statistics I: quantitative genetic parameters 

Common statistical analyses were performed using Statis- 
tica 6.0 software (2006). Our experimental design based on 
clonal genotypes, with the trait of interest assayed in the 
progeny, allowed us to estimate two variance components: G, 
the (broad sense) genetic variance, and E, the specific envi- 
ronmental variance or residual variance. We used a restricted 
maximum-likelihood based method, the animal model pro- 
cedure (Boldman et al. 1995) implemented in the MTD- 
FREML software. Iterations were continued until the dif- 
ferences in successive likelihood were less than 0.00001. We 
iterated the analysis in the full model (GE) and several start- 



ing values were assayed to ensure the solution being actually 
a global maximum of the likelihood. Also, we tested the pure 
environmental model by constraining the G parameter of the 
full model to zero, which gives a new likelihood value. The 
statistical significance of the G component was assessed us- 
ing likelihood-ratio tests (LRTs), where critical x 2 is selected 
with degrees of freedom equal to the number of parameters 
dropped from the model (i.e., df = 1 ; x 2 exit = 3.84 1 ). Genetic 
variances and covariances were estimated in a similar fashion, 
but including pairs of traits in a multivariate version of the an- 
imal model software (Boldman et al. 1995). We also present 
asymptotic standard errors (i.e., estimated values for large 
samples). 

Statistics II: G-matrices comparison 

The experimental treatment was sexual/asexual genotypes. 
The resulting G matrices were first compared using the Flury 
hierarchical method, using the software CPCrand (Phillips 
and Arnold 1999). This procedure is a principal component 
approach to the comparison of matrices (Flury 1988), whose 
application to G matrices was developed by Phillips and 
Arnold (1999). The method, based on maximum likelihood, 
compares two or more matrices and determines their struc- 
tural differences by comparison of eigenvectors (orientation 
of axes of the principal component analysis) and eigenvalues 
(variance along each axis) in a hierarchical fashion. This pro- 
gression tests (following this order): (1) unrelated structure, 
indicating that matrices do not share any eigenvector, (2) par- 
tial common principal component, indicating that matrices 
share some eigenvectors, (3) common principal components, 
indicating that matrices share all eigenvectors, but not eigen- 
values, (4) proportionality, indicating that matrices share all 
eigenvectors and eigenvalues all differing by a scalar amount, 
and (5) equality, indicating that matrices share all eigenvector 
and eigenvalues (see also Steppan et al. 2002). An LRT was 
calculated for each model against the model of "unrelated 
structure" (the "jump-up" approach; see Phillips and Arnold 
1999). A randomization was performed to test the null hy- 
pothesis that the model fits the data significantly better than 
the unrelated structure. We performed 5000 iterations, each 
randomly assigning genotypes to the compared group. We 
also compared the G matrices by the Jackknife-MANOVA 
method, which was performed by a script provided by 
D. Roff, and run in the S-plus software (version 2000). This 
procedure, developed by Roff (Roff 2002a, 2006), uses the 
Jackknife procedure to produce a distribution of pseudoval- 
ues of matrix elements within each group. These pseudoval- 
ues are produced by deleting each sampling unit (genotypes, 
in this case) in turn. The final data matrix is arranged such 
that the columns comprise the pseudovalues of each covari- 
ance and the rows are the results of the deletion of a given 
genotype (Roff 2006). 
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Results 

Trait means 

Pairwise comparisons between sexual and asexual genotypes 
showed that production of apterous females, lifetime repro- 
duction, and the intrinsic rate of increase (r) were signif- 
icantly higher in asexual compared with sexual genotypes 
(i.e., reproductive specialization; see Table 1). Both types of 
reproductive variants showed important inter- and intra- 
clone variation in all measured traits (Figs. 2 and 3). 

Heritabilities and genetic correlations 

In general, we found low levels of genetic variances (heri- 
tabilities being less than 0.20) in both reproductive variants, 
being somewhat higher in the asexual genotypes, where only 
SMR and the intrinsic rate of increase (r) resulted significant 
(Table 2). No significant heritability was found in sexuals 
and only three were found in asexuals: SMR, body size, and 
r. The fact that a genetic correlation existed between SMR 
and body size (see below) suggests that the heritability in 
SMR is a consequence of such correlation. Indeed, comput- 
ing heritability using body mass as a covariable yields low and 
nonsignificant heritability (data not shown). In the pooled 
sample, age at maturity also showed significant heritability 
(Table 2). 

The genetic correlations that resulted significant either for 
the whole sample or sexual and/or asexual genotypes are 
graphed in Figs. 2 and 3 (see also Table 3). We found seven 
large and positive genetic correlations in the asexual geno- 
types, which involved SMR, Mb, lifetime reproduction, apter- 
ous individual production, and r (Table 3). This contrasted 
with the outcome of sexual genotypes, where large intraclone 
variation produced only nonsignificant genetic correlations 
(Table 3; Figs. 2 and 3). Among them, the large positive 
genetic correlation between Mb and both fitness measures 
in asexuals (Fig. 2) suggests a strong association between 
body size and reproductive output in asexual genotypes. Also, 



both fitness measures showed a positive genetic correlation in 
asexual genotypes (Fig. 3) (Table 3). 

Matrix comparisons 

The differences between G matrices were not large, as we 
did not find any change in the sign of genetic covariances. 
Whereas the Flury method of P-matrix comparison yielded 
nonsignificant differences, the G-matrix comparison found 
that matrices shared their principal component structure 
(Table 4). However, the Jackniffe-MANOVA method indi- 
cated that matrices differed significantly (Table 4). Patterns 
of differentiation of asexual genotypes were explained by 
their larger values of genetic variance in SMR, fitness mea- 
sures and its genetic correlations with Mb, total fecundity, 
and apterous parthenogenetic females, compared with sex- 
ual genotypes. However, the P-matrix comparison did not 
show significant differences (Table 4). 

Discussion 

A general indicator of how a population is adapting to the en- 
vironment (assuming that the environment does not change) 
derives from the FTNS: the rate of increase in mean fitness 
in a population equals its genetic variance in fitness (Fisher 
1930; Gustafsson 1986; Frank and Slatkin 1992). In other 
words, the genetic variance in fitness of a given population 
would represent an instantaneous indicator of its mean rate 
of increase. In our experiment, we compared the quantita- 
tive genetics of life-history traits of sexual and asexual geno- 
types of an aphid in environmental conditions that promote 
sustained parthenogenesis. According to the FTNS, and as- 
suming that evolutionary equilibrium is not achieved (i.e., 
the population is always below the optimum), we predicted 
that the rate of increase in mean fitness in asexual genotypes 
should be higher than in sexual, which would be translated 
in heritability of fitness being significantly different than zero 
in asexuals, and higher than in the sexually specialized geno- 
types. Our results support these predictions (when using the 



Table 1. Phenotypic means (±SE from genotype means) of traits measured in asexual and sexual genotypes of R. padi in conditions of sustained 
parthenogenetic reproduction (see Methods for details). Nested-ANOVA was performed on clonal lines nested into each kind of reproductive mode. 
Significant effects were highlighted in bold and indicated by: *P < 0.05; "P < 0.01 ; ">P < 0.001 . 



Trait 


Asexual genotype 


Sexual genotype 


Clonal line 


Sexual/asexual 




Mean (±SE) 


Mean (±SE) 


F (33,242) 


F (1,242) 


Standard metabolic rate, SMR (ml_ C0 2 Ir 1 ) 


0.25 ± 0.009 


0.22 ±0.007 


2.044" 


1.925 


Adult body mass, M b (mg) 


0.59 ±0.011 


0.57 ±0.010 


1.947" 


1.531 


Age at maturity, AM (days) 


8.33 ± 0.078 


8.45 ±0.074 


1.947" 


0.350 


Apterous females, AF (individuals) 


54.0 ± 1.237 


49.7 ±1.269 


1.555* 


6.305* 


Winged females, VL (individuals) 


1.95 ± 0.434 


1.47 ±0.390 


1.321 


0.716 


Lifetime reproduction, LR (individuals) 


55. 9± 1.170 


51.2 ±1.232 


1.403 


7.190" 


Intrinsic rate of increase, r 


0.27 ± 0.003 


0.25 ± 0.003 


2.236 " 


18.22" 
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0.5 
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0.8 



=0.9910.29 




B 



0.4 



0 5 



0 6 



0.7 



0.8 



• Asexual 

=0.77±0.24 




0.4 0.5 0.6 0.7 0.8 

body mass 

Figure 2. Clonal means in sexual and asexual genotypes of ft. padi for: 
body mass and (A) standard metabolic rate (SMR), (B) total fecundity, 
(C) fitness (per-capita intrinsic growth rate). The genetic correlations for 
sexual and asexual genotypes are presented as a reference (see detailed 
statistics including the genetic correlation for the whole population in 
Table 3). 

intrinsic rate of increase as a measure of fitness), suggesting 
that asexual genotypes are being adapted to environmental 
conditions that promote parthenogenetic reproduction, rel- 
atively to the sexual genotypes. Similar results were found in 



0.36 
0.34 
0.32 
0.30 
0.28 
0.26 
0.24 
0.22 
0.20 



0.34 
0.32 
0.30 
0.28 
0.26 
0.24 
0.22 
0.20 




40 50 60 70 

apterous females 



• Asexual 
O Sexual 



r =0.98±0.19 




B 



40 50 60 70 

lifetime reproduction 

Figure 3. Clonal means in sexual and asexual genotypes of ft. padi for: 
(A) apterous female production and fitness; and (B) total fecundity and 
fitness. The genetic correlations for sexual and asexual genotypes are 
presented as a reference (see detailed statistics including the genetic 
correlation for the whole population in Table 3). 



rotifers (Becks and Agrawal 20 1 1 ), for which a trend seems to 
exist for a short-term disadvantage in facultative sexual lin- 
eages relative to their asexual counterparts. In more general 
terms, supports for the predictions of the FTNS come from 
plants, where locally adapted populations exhibited higher 
heritabilities of fitness, compared with those from the edge 
of the range (Emery et al. 20 1 1 ) , and also from evolutionary 
experiments in bacteria, which showed an increase in genetic 
variation in fitness in time, as adaptation to local environ- 
ments proceeds (Lenski and Travisano 1994). 

The structure of genetic covariances showed that asex- 
ual genotypes of R. padi exhibited several positive genetic 
correlations among traits, which were absent in the sexual 
genotypes. This result adds to a growing record of absence 
of life-history trade-offs when the contrary is predicted by 
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Table 2. Broad-sense heritabilities of a suite of life-history traits computed on sexual and asexual genotypes of ft. padi and measured in conditions of 
sustained parthenogenetic reproduction (see Methods for details). Significant values after REML and likelihood-ratio test with one degree of freedom 
are indicated with asterisks (/ 2 values are presented in parenthesis with: *P < 0.05, "P < 0.01, and "*P < 0.001) and highlighted in bold. 



Trait 



Asexual genotypes (W : 
REML (±SE) 



15) 



Sexual genotypes (W : 
REML (±SE) 



20) 



Pooled genotypes (W = 35) 
REML (±SE) 



Standard metabolic rate (SMR) 
Adult body mass (M b ) 
Age at maturity (AM) 
Apterous females (AF) 
Winged females (VL) 
Lifetime reproduction (LR) 
Intrinsic rate of increase (r) 



0.18 ±0.08 (10)" 
0.14 ± 0.08 (5.29)* 

0.08 ± 0.06 (1.74) 
0.06 ± 0.06 (3.33) 
0.03 ± 0.04 (0.90) 
0.06 ±0.06 (1.31) 
0.19 ±0.09 (3.91)* 



0.001 ±0.05 (0.01) 
0.06 ±0.06 (0.23) 
0.11 ±0.073 (0.09) 
0.04 ±0.05 (0.07) 
0.04 ±0.05 (0.89) 
0.01 ±0.05(1.12) 
0.08 ± 0.06 (0.87) 



0.11 ±0.05(7.83)** 
0.10 ±0.05(7.12)** 
0.10 ±0.05(6.23)* 

0.07 ±0.04 (3.44) 
0.04 ±0.04 (1.07) 
0.07 ±0.04 (3.62) 
0.17 ±0.05(16.21)** 



Table 3. Genetic correlations computed with REML on a suite of traits measured in sexual and asexual genotypes of ft. padi in conditions of sustained 
parthenogenetic reproduction (see Methods for details). Significant values after REML and likelihood-ratio test with one degree of freedom are 
indicated with asterisks (x 2 values are presented in parenthesis with: *P < 0.05, *"P < 0.01, and ***P < 0.001) and highlighted in bold. Abbreviations 
are as in Tables 1 and 2. 



Trait 


Asexual genotypes 


Sexual genotypes 


Pooled genotypes 




(±SE) 


REML (±SE) 


REML (±SE) 


SMR-M b 


0.99 ±0.11 (12.47)** 


0.01 ± 0.1 (0.23) 


0.99 ±0.11 (12.61)*** 


S MR-AM 


-0.39 ± 0.47 (0.59) 


0.12 ± 0.1 (1.46) 


-0.43 ±0.32 (1.49) 


SMR-APT 


0.19±0.09 (3.11) 


0.02 ± 0.2 (0.08) 


0.87 ± 0.35 (4.74)* 


SMR-VL 


-0.05 ±0.68 (0.08) 


0.04 ±0.1 (0.57) 


0.07 ±0.51 (0.02) 


SMR-LR 


0.94 ± 0.37 (3.86)* 


0.001 ±0.1 (0.10) 


0.91 ± 0.35 (5.25)* 


SMR-r 


0.66 ± 0.28 (3.91)* 


0.05 ±0.6 (0.14) 


0.70 ±0.21(5.79)* 


M b -AM 


-0.10 ±0.51 (0.03) 


-0.55 ±0.48 (0.01) 


-0.42 ±0.32 (1.39) 


M b -APT 


0.85 ± 0.27 (3.12) 


0.86 ± 0.68 (0.03) 


0.98 ± 0.21 (6.93)* 


M b -VL 


0.07 ±0.77 (0.11) 


0.66 ± 0.84 (0.35) 


0.66 ± 0.59 (1.55) 


M b -LR 


0.99 ± 0.29 (4.01)* 


1.0± 1.79 (0.11) 


1.0 ±0.31 (0.03) 


M b -r 


0.77 ± 0.24 (3.87)* 


0.85 ± 0.40 (0.05) 


0.81 ±0.17(8.10)** 


AM-APT 


-0.10 ± 0.64 (0.23) 


-0.81 ± 0.66 (1.98) 


-0.60 ±0.39 (2.01) 


AM-VL 


-0.81 ± 0.65 (1.47) 


-0.61 ± 0.64(0.83) 


-0.97 ± 0.53 (4.35)* 


AM-LR 


0.40 ± 0.60 (0.37) 


-1.0 ± 1.26 (0.03) 


-0.88 ± 0.36 (4.66)* 


AM-r 


-0.70 ± 0.29 (2.48) 


-0.44 ±0.44 (0.19) 


-0.57 ± 0.22 (3.80)* 


APT-VL 


-0.42 ± 0.80 (0.23) 


1.0 ±0.62 (2.11) 


0.01 ±0.62 (0.18) 


APT-LR 


0.97 ± 0.05 (1.63) 


1.0 ±0.35 (1.41) 


0.01 ± 1.56(0.22) 


APT-r 


0.84 ± 0.22 (3.92)* 


0.75 ± 0.39 (0.89) 


0.93 ±0.16 (8.07)** 


VL-LR 


-0.20 ± 0.95 (0.04) 


-1.0± 1.17 (0.01) 


0.001 ±0.60 (0.01) 


VL-r 


0.25 ± 0.67 (0.14) 


0.60 ± 0.80 (0.66) 


0.66 ±0.49 (2.21) 


LR-r 


0.98 ±0.19 (4.56)* 


1.0 ± 0.93 (0.78) 


1 .0 ± 0.99 (0.99) 



life-history theory (Vorburger 2005; von Burg et al. 2008; 
Castaneda et al. 2010c; King et al. 2011). In general, these 
evidences come from populations that are locally adapted 
or ecologically specialized, which exhibit genetic variation in 
resource acquisition. Hence, it appears that in addition of 
having comparatively high genetic variation in fitness, pos- 
itive genetic correlations among life-history traits are a key 
feature of locally adapted populations. Interestingly, these dif- 
ferences in genetic (co)variances between sexual and asexual 
genotypes were not evident after comparing the phenotypic 
(co)variance matrices but after the G-matrix comparison, 
which indicated that one principal component was shared. 



This is probably due to the fact that there were no differences 
in signs in the genetic correlations of sexual and asexual geno- 
types. 

As a measure of fitness, we used two common metrics 
for fitness determination in invertebrates: the intrinsic rate 
of increase (r) and lifetime reproduction (Huey and Berri- 
gan 2001). Both measures are supposed to be appropriate 
if the population is at mutation-selection equilibrium (Roff 
2010). However, when a population is in a growing phase, 
which is the case in our environmental conditions of sus- 
tained parthenogenesis, Roff (2010; see also Roff 2002b) has 
suggested that the best measure of fitness is r. In our dataset, 
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Table 4. G-P matrices comparisons between asexual and sexual lineages of R. padi using the Flury hierarchy and the MANOVA methods. The verdict 
of the Flury hierarchy is the model that best explains the difference between two matrices. The probabilities given for the Flury hierarchy correspond 
to the test of equality of two matrices, not necessarily to the verdict. The method MANOVA refers to the results obtained by using the scores of 
the pseudovalues of the G and P matrices. The numbers in parentheses for the MANOVA method correspond to the freedom degree associated to 
respective test. 





Flury 
Hierarchy 






MANOVA 
Jackknife 






LRT 


P-value 


Verdict 


X (Wills) 


F-value 


P-value 


G matrix 














Asexual-sexual 


3.116 


<0.05 


PCPC 2 


0.045 


17.06, 2 ,,, 7) 


<0.001 


P matrix 














Asexual-sexual 


5.647 


0.80 


Equal 


0.905 


1 -2/(21,255) 


0.19 



both measurements were strongly and genetically correlated 
only in asexual genotypes. Also, only r showed significant 
heritability in asexual genotypes (our interpretation as a con- 
sequence of the FTNS). It is reasonable to conclude, hence, 
that lifetime reproduction is perhaps not the best surrogate of 
fitness in experimental studies involving asexual reproducing 
organisms. 

Few studies have evaluated SMR variation in hemipterans 
(Salvucci and Crafts-Brandner 2000; Castaneda et al. 2010a,b; 
Artacho et al. 2011), which in general have not found spe- 
cific patterns associated with life-history strategies. Studies in 
other arthropods, however, suggest that increased reproduc- 
tive output could be associated with higher energy expendi- 
ture because of the maintenance of the metabolic machinery 
for reproduction (Terblanche et al. 2004; Lardies and Bozi- 
novic 2008). When analyzed as per gram basis, we did not 
find differences in energy expenditure between sexual and 
asexual genotypes. However, both energy expenditure and 
body mass showed positive genetic correlations with fitness 
traits, which suggests that asexual genotypes invest all their 
energy in the asexual reproduction, contrary to sexual geno- 
types that would be investing most of their energy in the 
sexual episode. 

According to some authors (Lee 2002; Bacigalupe 2009), 
the invasion success of many species might depend more 
heavily on their ability to respond to natural selection than 
on plasticity in physiological tolerance. In other words, it 
would be expected that the genetic variance of life histo- 
ries in invasive species is higher than in noninvasive species. 
This is not what the many invasive species that reproduce 
clonally suggests (Thompson and Eckert 2004; Voisin et al. 
2005; Kim et al. 2010). Several species of aphids have lost ge- 
netic variation in populations within the invaded range (e.g., 
the Americas; see Castaneda et al. 2010a,b) mainly because 
of the reduction in the capacity for sexual reproduction. In 
other words, specialized genotypes for asexual reproduction 
are the most successfully invasive genotypes as they maxi- 
mize fecundity at the expenses of survival insurance during 



winter. Here, the temperate climate is a key element in the 
invasive capacity of asexual genotypes. Our results confirm 
these observations, indicating that specialized asexual geno- 
types maximize fitness by changing their energy allocation 
scheme. 
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Table A1. Details of the 35 sexual and asexual genotypes of Rhopalosiphum padi used in this study, including the site, the date, and the plant of 
collection. 



Locus 



Clone code Collection site 



Collection date Plant 



R10.172 R5.10 R5.29b R5.50 



R6.27 



R3.171 S16b 



S17b 



asexl 


Pleine-Fougeres 


12/10/2007 


Wheat 


195,199 


256,260 


asex2 


Pleine-Fougeres 


28/09/2007 


Wheat 


195,199 


256,271 


asex3 


Pleine-Fougeres 


28/09/2007 


Wheat 


195,201 


258,260 


asex4 


Pleine-Fougeres 


12/10/2007 


Wheat 


195,201 


271,260 


asex5 


Pleine-Fougeres 


28/09/2007 


Wheat 


201,195 


254,264 


asex6 


Pleine-Fougeres 


28/09/2007 


Wheat 


197,201 


260,260 


asex7 


Pleine-Fougeres 


28/09/2007 


Wheat 


195,199 


258,260 


asex8 


Pleine-Fougeres 


28/09/2007 


Wheat 


201,201 


262,271 


asex9 


Pleine-Fougeres 


12/10/2007 


Wheat 


195,201 


271,260 


asexl 0 


Pleine-Fougeres 


28/09/2007 


Wheat 


195,197 


271,271 


asexl 1 


Pleine-Fougeres 


12/10/2007 


Wheat 


195,201 


271,260 


asexl 2 


Pleine-Fougeres 


20/1 1/2006 


Wheat 


195,201 


258,256 


asexl 3 


Pleine-Fougeres 


20/1 1/2006 


Wheat 


195,203 


256,260 


asexl 4 


Pleine-Fougeres 


20/1 1/2006 


Barley 


195,201 


271,256 


asexl 5 


Pleine-Fougeres 


20/1 1/2006 


Barley 


195,199 


256,271 


sexl 


La Chapelle Thouarault 


14/04/2008 


P. padus 


195,201 


262,256 


sex2 


La Chapelle Thouarault 


14/04/2008 


P. padus 


199,201 


260,271 


sex3 


La Chapelle Thouarault 


14/04/2008 


R padus 


195,201 


271,271 


sex4 


La Chapelle Thouarault 


14/04/2008 


R padus 


199,201 


260,271 


sex5 


La Chapelle Thouarault 


14/04/2008 


P. padus 


199,201 


254,271 


sex6 


La Chapelle Thouarault 


14/04/2008 


P. padus 


195,201 


254,256 


sex7 


La Chapelle Thouarault 


14/04/2008 


P. padus 


199,201 


266,271 


sex8 


La Chapelle Thouarault 


14/04/2008 


P. padus 


201,201 


264,271 


sex9 


La Chapelle Thouarault 


14/04/2008 


P. padus 


199,201 


262,271 


sex10 


La Chapelle Thouarault 


30/04/2008 


P. padus 


201,201 


260,262 


sexl 1 


La Chapelle Thouarault 


30/04/2008 


P. padus 


201,201 


260,264 


sex 12 


La Chapelle Thouarault 


30/04/2008 


P. padus 


199,201 


260,260 


sex 13 


La Chapelle Thouarault 


30/04/2008 


P. padus 


199,201 


271,260 


sex 14 


La Chapelle Thouarault 


30/04/2008 


P. padus 


199,195 


256,258 


sexl 5 


La Chapelle Thouarault 


30/04/2008 


P. padus 


201,195 


258,256 


sex 16 


La Chapelle Thouarault 


05/05/2008 


P. padus 


199,201 


271,271 


sex 17 


La Chapelle Thouarault 


05/05/2008 


P. padus 


199,201 


258,260 


sex18 


La Chapelle Thouarault 


05/05/2008 


P. padus 


201,201 


264,271 


sex19 


La Chapelle Thouarault 


05/05/2008 


P. padus 


195,201 


262,256 


sex20 


La Chapelle Thouarault 


05/05/2008 


P. padus 


199,201 


258,271 



172,181 
172,181 
178,179 
172,181 
174,185 
174,178 
174,174 
181,200 
174,174 
172,181 
172,202 
172,185 



172, 
172, 
195, 
172, 



172,204 
164,166 
181,181 
181,198 
172,185 
174,181 
181,185 
174,181 
166,172 
174,181 
172,200 
172,181 
172,208 
172,166 
172,181 
164,174 
181,174 
181,181 
178,181 



304,304 
304,328 
312,342 
304,336 
304,338 
324,340 
328,336 
304,310 
304,312 
304,366 
304,314 
304,304 
304,360 
304,314 
304,328 
304,328 
304,336 
304,330 
304,310 
304,310 
304,304 
304,322 
304,310 
328,338 
308,347 
304,338 
304,342 
304,304 
304,334 
304,356 
304,310 
304,312 
304,338 
304,304 
304,350 



178,195 
178,195 
188,188 
178,195 
178,184 
180,184 
190,191 
178,180 
178,186 
178,195 
186,195 
195,195 
178,195 
190,195 
178,195 
195,211 
190,195 
178,178 
178,180 
178,178 
184,195 
178,178 
180,184 
178,178 
186,195 
178,178 
178,195 
180,195 
178,195 
195,195 
180,195 
178,180 
178,180 
178,180 
178,193 



225,227 
216,225 
218,222 
214,225 
207,225 
214,225 
216,216 
225,247 
216,225 
225,225 
225,243 
225,225 
225,225 
227,243 
225,225 
227,229 
216,225 
216,225 
227,247 
225,247 
207,225 
225,243 
225,247 
207,216 
216,243 
207,225 
227,227 
207,227 
216,225 
216,227 
225,247 
207,225 
225,225 
225,227 
225,235 



147,149 
147,149 
147,147 
147,149 
147,147 
147,147 
147,147 
147,147 
147,147 
147,149 
147,149 
147,149 
147,149 
147,149 
147,149 
147,149 
147,147 
147,147 
147,147 
147,147 
147,149 
147,147 
147,147 
147,147 
147,147 
147,147 
147,147 
147,147 
147,149 
147,149 
147,147 
147,147 
147,147 
147,149 
147,147 



159 
161 
159 
159 
159 
161 
159 
159 
159 
159 
159 
159 
165 
159 
159 
159 
159 
159 
159 
159 
161 
165 
171 
161 
159 
159 
159 
159 
159 
159 
165 
159 
159 
159 
159 
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Figure A1. Neighbour-joining tree based on the allele shared distance (DAS, Chakraborty and Jin 1993) matrix calculated with eight microsatellite 
loci for the 35 sexual and asexual multilocus genotypes of Rhopalosiphum padi used for fitness and energetic measurements. 
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